Astronomy & Astrophysics manuscript no. 


February 1, 2008 


(DOI: will be inserted by hand later) 





Near- 1 R spectroscopy of planetary nebulae precursors * 

D.A. Garcia-Hernandez 1 , A. Manchado 1 ' 2 , P. Garcfa-Lario 3 , C. Dommguez- Tagle 1 , G.M. Conway 4 , and F. 

Prada 5 

1 Instituto de Astrofi'sica dc Canarias, E-38205 La Laguna, Tenerife, Spain 
e-mail: (agarcia, amt , carlosd)@ll . iac . es 
Consejo Superior de Investigaciones Cientfficas, CSIC 

3 ISO Data Centre, Science Operations and Data Systems Division, Research and Scientific Support Department 
of ESA, VILSPA, E-28080 Madrid, Spain 

e-mail: pgarcia@iso . vilspa.esa.es 

4 Department of Physics and Astronomy. University of Calgary, Alberta T2N 1N4, Canada 

5 Centro Astronomico Hispano— Aleman. Apartado de Correos 511, E-04080, Almerfa, Spain 

Received July 31, 2001; accepted March 18, 2002 

Abstract. We present near-IR spectroscopy of a sample of 30 IRAS sources recently identified as late AGB stars, 
post-AGB stars or early PNe. The spectra obtained are centered at various wavelengths covering the molecular 
hydrogen v=l— »Q S(l) 2.122 fim and v=2— >1 S(l) 2.248 /jm emission lines, the recombination lines of hydrogen 
Br7 2.166 /im, Pfy 3.741 /im and Bra 4.052 fim, and the CO[i>=2— >0] first overtone bandhead at 2.294 /jm. 
As a result of these observations we have detected molecular hydrogen emission for the first time in 9 of these 
sources and confirmed a previous detection by Weintraub et al. (1998). This increases from 4 to 13 the total 
number of proto-PNe detected in molecular hydrogen. In most cases, the positive detections also show emission 
in the recombination lines of hydrogen (with the exception of IRAS 17150—3224) indicating that the onset of 
molecular hydrogen emission takes place in the post-AGB phase, very shortly before the nebula becomes ionized. 
When the molecular hydrogen is fluorescence-excited the detection rate is found to be directly correlated with 
the evolutionary stage of the central star, rather than with the nebular morphology. When the temperature of 
the central star is hot enough, fluorescence excitation can be induced by the absorption of UV photons escaping 
from the rapidly evolving central post-AGB star. In contrast, shocked-excited molecular hydrogen is detected only 
in strongly bipolar proto-PNe, sometimes even at an early stage in the post-AGB phase. Shock-excitation is the 
consequence of the interaction of the fast post-AGB wind with the slow wind material ejected during the AGB. 
The strong correlation of shocked-excited molecular hydrogen emission with bipolarity found confirms the result 
previously reported by Kastner et al. (1996) in evolved PNe. However, our results show that this correlation does 
not exist in the case of fluorescence-excited molecular hydrogen. 

Key words. Stars : AGB and post-AGB - planetary nebulae: general - stars: evolution - Infrared: stars - stars: 
circumstellar matter - stars: mass loss 



1. Introduction 

The short transition phase between the end of the asymp- 
totic giant branch (AGB) and the formation of a new plan- 
etary nebula (PN) is still poorly understood. The main 
reason for this is probably twofold. First, only a small 
number of objects are known in this phase, as this is 
a short-lived evolutionary stage. Second, in many cases, 
stars evolve through this phase heavily obscured by the 
strong circumstellar envelopes formed during the previous 
AGB phase, thus making observations difficult. 
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However, this is an evolutionary stage when crucial 
changes occur that completely determine the subsequent 
evolution of the individual star as a PN. Eventually, this 
also results in changes in the overall chemical evolution of 
the Galaxy, since huge amounts of processed material are 
transferred to the interstellar medium as a consequence of 
the strong mass loss experienced at the very end of the 
AGB phase. 

It is also during this short transition phase when the 
ionization of the gas ejected during the previous AGB 
phase takes place. However, it is now well established that 
in most cases only a small fraction of it gets eventually 
ionized in comparison with the total amount of mass ex- 
pelled by the central star during the AGB phase (tenths 
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of a solar mass of ionized gas compared to up to several 
solar masses ejected during the AGB phase in the case of 
the more massive progenitors) (Pottasch 1980, 1992; Kwok 
1994). Most of this material remains neutral in the form 
of dust grains, molecules or atoms which can be detected 
easily at infrared or radio wavelengths even when the star 
has completely evolved and turned into a PN. 

From the analysis of IRAS data, and especially with 
the advent of the higher quality ISO data, it has been 
possible to confirm the presence of large amounts of dust 
grains which emit thermally in the mid- and far-infrared in 
objects evolving from the AGB to the PN stage. Through 
the analysis of the features observed in the ISO spectra it 
is possible to determine the dominant gas and dust chem- 
istry (O-rich or C-rich) in the stellar photosphere and in 
the circumstellar envelope (Waters et al. 1999; Sylvester 
et al. 1999; Hrivnak et al. 2000). 

The detection of molecules in the envelopes around 
PNe has also been successful at submilimeter and radio 
wavelengths over the last two decades. Strong CO molecu- 
lar emission has been detected in a large number of C-rich 
PNe (Huggins & Healy 1989; Huggins et al. 1996) and 
is also detectable around both C-rich and O-rich post- 
AGB stars (Likkcl et al. 1991; Loup et al. 1993; Neri et al. 
1998, and references therein). Other molecules, like HCN 
or CN are detected only in envelopes dominated by a C- 
rich chemistry (Lindqvist et al. 2000), while in the case 
of circumstellar envelopes dominated by an O-rich chem- 
istry strong molecular maser emission from SiO, H2O and 
OH is usually detected when the star is in the AGB phase 
(Engels 1979; Engels & Lewis 1996; Nyman et al. 1998; te 
Lintel Hekkert 1991a; te Lintel Hekkert et al. 1991b). In 
many cases, the OH maser emission remains detectable be- 
yond the AGB throughout the whole post- AGB evolution 
(te Lintel Hekkert & Chapman 1996) and sometimes even 
when the object has already become an ionized PN (the so- 
called OHPNe - Zijlstra et al. 1989; 1991). More complex 
molecules are difficult to detect, but in the circumstellar 
envelopes of a few objects like IRC+10216, AFGL 618 
and AFGL 2688 more than 50 different molecular species 
have been identified in the last few years, and the number 
and complexity of the new molecules found is continuously 
increasing (Wallcrstein & Knapp 1998; Bujarrabal et al. 
1988; Cernicharo et al. 1999). 

Since hydrogen is the main constituent of stellar atmo- 
spheres one would expect to find also a significant amount 
of neutral material locked up in the form of molecular hy- 
drogen around stars evolving from the AGB phase to the 
PN stage (both O-rich and C-rich). Following the first 
detection of molecular hydrogen in the PN NGC 7027 
(Treffers et al. 1976), more than 70 other PNe have been 
detected in the last few years by Webster et al. (1988), 
Zuckerman & Gatley (1988), Kastner et al. (1996), Latter 
et al. (1995), Hora et al. (1999) and Guerrero et al. (2000). 
As a general trend, it is found that molecular hydrogen 
emission is predominantly detected in type I bipolar PNe, 
which are expected to be the result of the evolution of the 



more massive AGB progenitors. It remains unclear, how- 
ever, whether this is just an observational selection effect. 

Very few attempts have been made so far to search for 
molecular hydrogen around post- AGB stars or proto-PNe, 
the precursors of PNe. Molecular hydrogen emission was 
first detected towards AFGL 618 by Beckwith et al. (1978) 
and AFGL 2688 by Gatley et al. (1988). More recently, 
Weintraub et al. (1998) looked at a small sample of five 
proto-PNe and detected molecular hydrogen emission in 
two of them (IRAS 17441-2411 and AFGL 6815S=IRAS 
17150-3224), bringing to four the total number of bipo- 
lar proto-PNe with this type of emission. The discovery 
and study of new transition objects showing this emis- 
sion at such an early stage is crucial for determining the 
mechanism of excitation of molecular hydrogen and the 
conditions needed to preserve this neutral material from 
destruction in more evolved PNe. 

This paper extends this search to a larger sample of 
transition objects, including late-AGB stars, post-AGB 
stars with various spectral types, as well as the youngest 
PNe known, in order to look for possible correlations be- 
tween the onset of molecular hydrogen emission and the 
evolutionary stage and/or the morphology of the sources 
observed. 

In Sect. 2 we explain the criteria followed to select 
the sources in the sample. The spectroscopic observations 
made in the near infrared are described in Sect. 3, while 
the results obtained are presented in Sect. 4 and discussed 
in Sect. 5. The conclusions derived from these results are 
given in Sect. 6. 

2. Selection of the sample 

For the last ten years our group studied a large num- 
ber of IRAS sources with dust temperatures similar to 
post-AGB stars and planetary nebulae (^100—200 K), 
mainly through near-infrared photometry (Garcfa-Lario 
et al. 1990; 1997a; Manchado et al. 1989a) and optical 
spectroscopy (Garcfa-Lario et al. 1991; 1994; 1997b; 2001; 
Ricra et al. 1995). As a result of this multi- wavelength 
program we have discovered heavily obscured late AGB 
stars and early post-AGB stars still surrounded by thick 
circumstellar envelopes only detectable in the near in- 
frared as well as post-AGB stars with optically bright 
central stars covering all possible spectral types from M 
to B in what it seems to be an evolutionary sequence to- 
wards the PN stage, together with a small number of new 
extremely young PNe (Garcfa-Lario et al. 1997a; 2002), 
among them, the youngest PN known so far: the Stingray 
Nebula (Bobrowsky et al. 1998). 

From this list of newly discovered sources, we selected 
a subsample of 30, which are observable from the south- 
ern hemisphere. This subsample is also representative of 
the various evolutionary stages mentioned above, includ- 
ing very late AGB stars with no optical counterparts, 
post-AGB stars with all kind of spectral types and proto- 
PNe showing a variety of morphological structures (from 
round/elliptical to strongly bipolar). 
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Table 1. IRAS sources included in the survey and main observational properties. 
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In Table 1 we list the IRAS sources included in the 
study. In this Table we also indicate some of the observa- 
tional characteristics which are relevant for the discussion 
presented in Sect. 5, like their apparent optical morphol- 
ogy, spectral type, or the presence/absence of Ha emission 
in their optical spectra (all this information taken from 
Garcfa-Lario et al. 2002, and references therein). 

3. Observations and data reduction 

The observations were conducted in March 1993 and 
January 1994 with IRSPEC attached to the Nasmyth fo- 
cus of the 3.5 m NTT telescope at La Silla, ESO. IRSPEC 
was a cryogenically cooled scanning grating IR spectrom- 
eter with a 58x62 pixel array from SBRC as detector 
(Moorwood et al. 1986). It covered the 1—5 microns wave- 
length range with a spectral resolution between 1500 and 
2000. The spatial scale was 2.2"/pixel, the slit width was 
4.5"and the un-vignetted field was ~100". The array had 
a quantum efficiency of 0.89 at 2.85 ^m, and a well capac- 
ity of 1 10 6 e~. This allowed exposures up to 60 seconds 
per frame in the 2—2.6 /xm region. For integration times 
around 10 seconds, the combined dark and read noise was 
^150 c~. Owing to the poor wavelength coverage of the 
spectrograph in a single exposure, different exposures were 
needed centered at various wavelengths, each one cover- 



ing a different spectral feature under analysis. The lines 
observed were the molecular ro- vibrational lines of molec- 
ular hydrogen u=1^0 S(l) 2.122 ^m and w=2^1 S(l) 
2.248 /im, the recombination lines of hydrogen Br7 2.166 
^m, Pf 7 3.741 fim and Bra 4.052 ^m and the CO[w=2^0] 
2.294 fim first overtone bandhead. Total on-source integra- 
tion times were typically of 240 s with exposures ranging 
from 36 s to 1440 s in the most extreme cases. 

The slit orientation was north-south (P.A. — 0°) in 
1993 with the exception of IRAS 10178-5958 (P.A. = 
70°), IRAS 16594-4656 (P.A. = 70°), IRAS 16342-3814 
(P.A. = -70°) and IRAS 17423-1755 (P.A. = +55°), 
while in 1994 the slit orientation was always east- west 
(P.A. = 90°). The changes in the slit orientation during 
the 1993 run were in principle made in order to follow the 
bipolar axis of some well known extended objects, like the 
ones listed above. However, due to a mistake in the de- 
termination of the sign, IRAS 17423—1755 was observed 
almost perpendicular to the actual bipolar axis in 1993 
while in IRAS 16342-3814 independently of the slit ori- 
entation, we cover the source completely. The only other 
known extended bipolar objects in our sample are IRAS 
10197-5750 and IRAS 17150-3224. IRAS 10197-5750 
was only observed in 1994, when the orientation of the 
slit was 90°, while the bipolar axis of this source is lo- 
cated at P.A. ~ 20°. This means that the slit was posi- 
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Table 2. Summary table with the results obtained. Quoted values are the emission line fluxes or 2a upper limits to 
the expected emission line fluxes (in units of 10~ 14 erg cm -2 s -1 ). The equivalent widths (in units of 1CP 4 /im) are 
shown in brackets for the features detected in absorption. 
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tioned also very close to the direction perpendicular to its 
bipolar axis. In the case of IRAS 17150—3224, we observed 
the source both in 1993 and 1994 with orientations of 
and 90 "respectively, while the actual bipolar axis of the 
source is now known to be located at P. A. ~ —45°. Two 
other sources (IRAS 10178-5958 and IRAS 16594-4656) 
also were observed twice. The implications of these mis- 
adjustments in the analysis of the results obtained will be 
discussed later. 

The standard beam— switching technique was used in 
order to subtract the sky background. The data reduc- 
tion process includes bias and flat-field corrections, as well 
as wavelength calibration, sky subtraction, and absolute 
flux calibration using reference standard stars. For this 



we made use of the S/W package IRSPEC installed on 
MIDAS (version November 1998). 

4. Results 

In Table 2 we list the detections/non-detections made as 
a result of the observations described in Sect. 3. 

All sources but one in the sample were observed in the 
wavelength region surrounding the H 2 [v=1^0 S(l)] line, 
resulting in positive detections in 10 cases. Out of them, 9 
sources were also observed in the weaker H 2 [w=2^1 S(l)] 
line, among which 5 showed some emission. 

The Br7 line was detected in 17 objects, either in emis- 
sion (11) or in absorption (6). In a few cases, observations 
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were also performed in other spectral lines like Pf 7 , and 
Bra, or around the CO[u=2^0] first overtone bandhead, 
in order to obtain complementary information. 

The quantitative results are also shown in Table 2. 
Quoted values are the emission line fluxes derived by fit- 
ting gaussians to the observed line profiles or 2a upper 
limits to the expected emission line fluxes in the case of 
those sources with no emission detected over the underly- 
ing near infrared continuum. Both are shown in units of 
10~ 14 erg cm~ 2 s _1 . We have also estimated the equivalent 
widths of those features detected in absorption (in units 
of 10 -4 lira). The values obtained are shown in Table 2 
within brackets. 

In a few cases, we did not detect any significant signal 
in the 2D spectra (nor even a continuum). These sources 
have been classified in Table 2 as not detected. In these 
latter cases we can generally infer the object's position 
on the 2D frames based on information obtained at other 
spectral lines and estimate the associated rms of the pixels 
where the feature investigated had been detected. Then we 
use these rms to calculate the quoted values as 2a upper 
limits to the expected emission line fluxes in Table 2. 

We estimate that for a typical 240 s total on-target 
exposure and when the underlying near infrared contin- 
uum is weak our 2a detection threshold is of the or- 
der of - 1(T 15 erg cm" 2 s- 1 for the U 2 [v=1^0 S(l)], 
Br 7 , R 2 [v=2^1 S(l)] and CO[w=2^0] lines and of - 
10 -14 erg cm -2 s -1 in the Pf 7 and Bra regions. Note that 
the 2a upper limits in Tabic 2 arc very different from 
source to source because the rms is dependent of the con- 
tinuum. 

Three objects (IRAS 10178-5958, IRAS 16594-4656 
and IRAS 17150-3224), were observed in both runs (1993 
and 1994) and for them we quote in Table 2 a mean value 
of the integrated line fluxes (the flux differences from 
epoch to epoch are less than 20%). However, the emis- 
sion detected was always found to be consistent with a 
point-like nature of the observed source with the excep- 
tion of IRAS 10178-5958 and IRAS 17119-5926, where 
the emission found is clearly extended (see Sect. 5.5). 

Spectra for all objects (except for IRAS 05284+1945, 
whose near infrared counterpart was not detected at any 
wavelength) in Br 7 , H 2 [v=1^0 S(l)], H 2 [v=2^1 S(l)], 
CO[u=2— >0], Pf7 and Bra are presented in Figures la to 
If. For those sources that were observed twice we show 
the best spectrum. 

5. Discussion 

5.1. Hydrogen recombination lines 

In general, the strength of the near infrared hydrogen lines 
detected in the stars in our sample is consistent with the 
spectral types derived from optical data. All objects where 
a Br7 absorption has been detected show Ha absorption 
in the optical or very weak emission, while those with Br 7 
strongly in emission in the near infrared are also strong 
Ha emitters in the optical, as expected. 



In the case of the late AGB stars included in our study 
(those classified as OH/IR stars in Table 1) the Br 7 line 
region is veiled by the strong dust continuum emission 
and we are not able to detect any signature of this hy- 
drogen recombination line. These stars are pulsating with 
very long periods so heavily obscured by their circumstel- 
lar envelopes that they usually do not show any optical 
counterpart. 

The same veiling effect might be the reason why IRAS 
08005—2356 shows a featureless continuum around the 
Br 7 line as well as at other near infrared hydrogen re- 
combination lines, while in the optical a relatively strong 
Ha emission is detected. Another possibility is that the 
emission is just strong enough to fill in the photospheric 
absorption and, thus, both components compensate each 
other, leading to a flat continuum. 

IRAS 09024—5019 seems to be also a remarkable case, 
since it shows a clear emission in the light of Br 7 and Bra 
while it does not have any (known) optical counterpart. 
This star could be a heavily obscured transition object in 
which the onset of the ionization is now taking place. 

5.2. CO first overtone 

The effective temperature of late AGB stars is expected 
to be very low (<2500 K). Consistent with this, a 
strong CO photospheric absorption was detected in IRAS 
08425-5116 at 2.29 /xm, the only late AGB star included 
in our sample for which we obtained a spectrum at this 
wavelength. 

This feature corresponds to the CO first overtone 
bandhead, which is expected to be visible in stars with 
effective temperatures below 5000 K. However, IRAS 
17150—3224, a post-AGB star with a cool central star 
(G2I), did not show any signature of this bandhead, prob- 
ably because of the veiling produced by the strong thermal 
emission from the hot dust formed in the circumstcllar en- 
velope, which completely dominates the observed spectral 
energy distribution at this wavelength. 

Noteworthy, in three of the post-AGB stars with the 
earlier spectral types in the sample we detected the CO 
first overtone bandhead in emission. This has also been ob- 
served in a few other transition objects before (Oudmaijer 
et al. 1995) and it has been interpreted as the result of ac- 
tive post-AGB mass loss, the CO lines being formed in 
the dense outflow. 

5.3. H2 emission and evolutionary stage 

In order to understand the mechanism of excitation of 
molecular hydrogen it is worth to try to find a possible 
correlation between the frequency of positive detections 
and the evolutionary stage (spectral type of the central 
star) of the sources observed. 

Starting from the less evolved objects, we find that 
none of the late-AGB stars included in the sample show 
any indication of the presence of H 2 emission at 2.122 /im. 
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Fig. la. Near infrared spectra of the sources detected with IRSPEC. Each row corresponds to observations of a different object. 
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It is important to note however that, as for the hydrogen 
recombination lines above discussed, the spectral region 
around this molecular hydrogen line appears usually com- 
pletely dominated by the thermal emission coming from 
the hot dust in the envelope, making the detection of any 
spectral feature difficult in this kind of objects. 

More interesting is the fact that, with the only excep- 
tion of IRAS 17150-3224, which will be discussed below, 
none of the post-AGB stars in the sample with spectral 
types later than A showed any detectable emission in this 
line. This includes 6 F-type, 2 G-typc, 1 K-type and 1 M- 
type post-AGB stars. It should be noticed that in these 
optically bright post-AGB stars the veiling effect from the 
circumstellar dust is small and cannot justify the non- 
detections. 

In contrast, all the stars identified as post-AGB can- 
didates in our list with spectral types earlier than A were 
detected in molecular hydrogen at 2.122 ^m. In this case, 
the only exception was IRAS 17423—1755, whose central 



star has tentatively been classified as A-type or B-type 
and it shows no evidence of H2 emission at 2.122 /zm. 
If this emission comes mainly from the bipolar lobes, it 
may have not been covered by our slit because it was 
positioned by mistake perpendicular to the bipolar axis 
(see Sect. 3). Consistently, we did not detect any extended 
emission from this object in any other spectral line (nor 
even in Br7 which is detected very strong in emission). 
In addition, this is also a quite peculiar, extended bipo- 
lar proto-PN which shows a very high velocity field and 
a strong near infrared excess (Riera et al. 1995; Garcfa- 
Lario et al. 1997a). Thus, the H 2 emission may be veiled 
by the strong underlying continuum and/or the molecu- 
lar hydrogen molecules may have been destroyed by the 
strong stellar winds. 

Out of the 10 objects in which we detected molecu- 
lar hydrogen emission, 9 are new detections (only IRAS 
17150—3224 was previously reported by Weintraub et 
al. 1998). This brings the total number of transition 
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stars in which molecular hydrogen emission has been de- 
tected from 4 to 13. Note that the H 2 emission in IRAS 
10178—5958 was previously presented by us (Garcfa-Lario 
et al. 1999b). 

The results above are consistent with the non- 
detection of molecular hydrogen emission in the bipo- 
lar proto-PNe OH 231.8+4.2, IRAS 07131-0147 or IRAS 
09371+1212 by Wcintraub et al. (1998), all of which have 
late-type M central stars. 

However, IRAS 17150-3224 is not the only post-AGB 
star with a spectral type later than A showing H 2 emis- 
sion in the short list of known H 2 emitters. There is an- 
other discrepant case in the literature, the well-known 
PPN CRL 2688 (Sahai et al. 1998), whose central star is 
classified as F5Iae, although it is suspected to be a binary. 

5.4. H2 emission and morphology 

It has been suggested that molecular hydrogen emission 
is a signpost for bipolarity in evolved PNe (Kastner et al. 



1996). In order to explore the validity of this statement 
also in the precursors of PNe we can search for a possible 
correlation between the detection rate of molecular hydro- 
gen emission and the morphology of the sources included 
in our sample. 

For this, we will restrict our analysis in the follow- 
ing to those sources for which morphological information 
exists. This means that we will not consider, of course, 
any of the heavily obscured objects in the list which do 
not show any optical counterpart (all the late-AGB stars 
plus IRAS 09024-5019), neither those objects classified 
as of unknown size/morphology in Table 1. Note that 
this class of objects may actually contain a mixture of 
round/elliptical and bipolar sources, as deduced from re- 
cent observations of similar objects made with HST (al- 
though they could look point-like sources in ground-based 
CCD images) (e.g. Ueta et al. 2000). 

The remaining sources in the sample were classified as 
round/elliptical or bipolar depending on the morpholog- 
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Fig. lc. Near infrared spectra of the sources detected with IRSPEC (cont.). The best ¥L2[v—l~ +0 S(l)] and Br7 spectra for 
IRAS 10178-5958 correspond to the 1994 run 



ical properties derived from ground-based or HST obser- 
vations (see again Table 1). 

If we concentrate our attention now on the ten sources 
showing H2 emission in our sample we can see that mor- 
phological information is available for six of them, out of 
which five are classified as bipolar in Table 1, while only 
one is considered to be round/elliptical. 

However, among the bipolar objects in our sample 
(eight in total), there are three of them which were not 
detected in H 2 . This argues against the idea of a one to 
one correlation between H 2 emission and bipolarity. 

In summary, our main conclusion is that although 
there seems to be an association between bipolarity and 
the presence of H2 emission the correlation is not so strong 
as previously reported by Kastncr et al. (1996), who de- 
rived this conclusion from the observations of well evolved 
PNe. 



5.5. The excitation mechanism of H2 

Molecular hydrogen can be excited either in shocks or by 
UV photons (Burton 1992). In order to investigate which 
of these two possible mechanisms is responsible for the 
emission observed in the objects included in our sample it 
is useful to compare the relative strength of the H 2 [v=l— >0 
S(l)] and H 2 [u=2^1 S(l)] lines. From the analysis of this 
diagnostic ratio it is possible to determine whether the 
excitation is produced by shocks (Burton, Hollcnbach & 
Tielens 1992) or by fluorescence through the absorption 
of UV photons coming from the central star (Black & van 
Dishoeck 1987; Sternberg & Dalgarno 1989). For moderate 
UV radiation fields and gas densities lower than 10 5 cm~ 3 , 
shock— excited emission is expected to be associated to 
flux ratios H 2 [v=1^0 S(l)/w=2^1 S(l)] - 10 while in 
the case of emission induced by fluorescence a value of ~ 
2 (Aspin et al. 1993; Ramsay 1993; Burton et al. 1998) 
is predicted. In Table 3 we show the ratios derived for 
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9 out of the 10 post-AGB stars in the sample detected 
in molecular hydrogen for which data is available. The 
estimated uncertainty is ^20%. For those objects in which 
the H 2 [v=2— >1 S(l)] line was not detected lower limits are 
given. 

As we can see in Table 3 a wide range of values is 
obtained. If the above mentioned models are applied to 
the stars in our list we deduce that shock— excited H 2 
emission must be the dominant excitation mechanism 
in at least three of the objects in our sample with a 
clearly marked bipolar morphology: IRAS 10178—5958 , 
IRAS 16594—46456 (both of uncertain spectral type) and 
IRAS 17150-3224 (G-type central star) as well as in two 
other objects of unknown morphology: IRAS 14331—6435 
and IRAS 17311-4924 (both with B-type central stars), 
while fluorescence excitation seems to be the preva- 
lent excitation mechanism in the Stingray Nebula (IRAS 
17119—5926), which is a very young PN showing an in- 
cipient bipolar morphology. For the remaining objects the 



lower limits found are ambiguous to support cither fluo- 
rescence or shock excitation. 

From the results obtained, there seems to be a preva- 
lence of shocked-excitation in those objects showing the 
most extreme bipolar morphologies. However, the results 
are not conclusive since the statistics are very poor. 
There is also no evident connection between the excitation 
mechanism derived from the analysis of the H 2 [u=1^0 
S(l)/w=2^1 S(l)] ratio and the spectral type of the cen- 
tral star in which molecular hydrogen has been detected. 
However, in this latter analysis we must be careful since 
it is well known that episodic mass loss during the post- 
AGB phase can induce significant spectral changes which 
may alter for some time the observable spectral type and, 
thus, taking the spectral type as an indicator of the evo- 
lutionary stage might be not the best choice (Bobrowsky 
ct al. 1998; Garcfa-Lario ct al. 2001). 

In addition, an important caveat needs to be consid- 
ered when interpreting the data in Table 3. This is the fact 
that the models above mentioned lose their validity at very 
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Fig. le. Near infrared spectra of the sources detected with IRSPEC (cont.). The best H2[u=1^0 S(l)] spectra for IRAS 
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high densities (n e > 10 5 cm 3 ) and we know that some H 2 [ 



of the post-AGB stars and transition objects under study 
are surrounded by regions of these very high densities. In 
this case, collisional de-excitation of fluorescence excited 
levels ('collisional fluorescence') might become important, 
increasing the v=0 and 1 populations over the pure fluo- 
rescent values, leading to higher H 2 [u=l— >0 S(l)/v=2— >1 
S(l)] ratios (Sternberg & Dalgarno 1989; Burton et al. 
1990) which could mimic shock-excitation. 

In order to address the problem from a different point 
of view it is also useful to study the detailed spatial distri- 
bution of the molecular hydrogen emission and compare 
this, whenever possible with the spatial distribution ob- 
served in the light of Br7. 

The H 2 [v=1^0 S(l)]/Br7 ratio has recently found to 
be correlated both with the evolutionary stage and the 
degree of bipolarity in a sample of well known bipolar 
PNe (Guerrero et al. 2000). Bipolar PNe with small size, 
poorly defined rings and bright central stars exhibit small 



H 2 [w=1^0 S(l)]/Br 7 ratios (~ 0.1-0.5) while those with 
large, well defined rings and faint central stars showed 
larger H 2 [w=1^0 S(l)]/Br7 ratios (>5). 

Guerrero et al. (2000) suggest that the excitation of the 
molecular hydrogen in each of these groups turns from flu- 
orescence excitation at the early stages of PNe to shocked- 
excited emission as the bipolar PNe evolve. This is sup- 
ported by the different spatial distribution of the observed 
molecular hydrogen emission. 

H2-dominated bipolar PNe are found to show different 
morphological and physical properties. Their central re- 
gions have low surface brightness, in agreement with their 
low densities while Br7-dominated bipolar PNe have a 
high surface brightness and larger densities. 

If this effect is interpreted as an evolutionary sequence 
(Balick 1987) then the less evolved bipolar PNe with 
molecular hydrogen emission would start their evolution 
as Br7 dominated becoming later H 2 -dominated. 
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Fig. If. Near infrared spectra of the sources detected with IRSPEC (cont.) 



Note, however, that this scenario docs not apply 
to all PNe. Well evolved, low-mass density-bounded 
round/elliptical PNe are not expected to preserve hydro- 
gen in the form of molecules since most (if not all) of the 
hydrogen in the circumstellar envelope becomes ionized as 
the ionization front passes through the gas. These PNe will 
end their lives as Br7-dominated, unlike the group of PNe 
with more massive progenitors. Those with more massive 
progenitors are expected to develop ionization-bounded 
nebulae, where a substantial fraction of the mass ejected 
will never become ionized. This enables the detection of 
neutral material in the circumstellar gas even at a very 
late stage as PNe, provided the UV radiation field from 
the central star is hard enough to produce fluorescence 
excitation. In the subgroup of bipolar PNe, the shielding 
provided by the dense equatorial regions that collimate 
the outflow from the UV photons escaping from the cen- 
tral star would also favor the preservation of hydrogen 
in molecular state. This would explain the high detec- 
tion rate of molecular hydrogen emission in bipolar type I 
PNe, which are suspected to be the result of the evolution 
of stars with large progenitor masses. In this latter case, 
molecular hydrogen can also be excited by shocks formed 
in the high velocity outflow. 



Applying the above considerations to the analysis of 
the molecular hydrogen emission observed in precursors of 
PNe it would be natural to expect fluorescence to be ac- 
tivated when the radiation coming from the central post- 
AGB star becomes hard enough to produce a significant 
number of UV photons as the central star moves in the 
HR diagram towards hotter effective temperatures in its 
way to becoming a PN. On the other hand, the interaction 
of the fast post-AGB wind with the slow-moving material 
ejected during the previous AGB phase would favor the 
formation of shocks, especially in those objects showing 
a strong collimation induced by the presence of a thick 
equatorial disk or torus, in the process of developing a 
bipolar morphology. 

In order to check whether this scenario is consistent 
with the data obtained on the stars in our sample we 
have derived the H2[u=1^0 S(l)]/Br7 ratios for the ob- 
jects that exhibit H2 emission. The results obtained are 
presented in Table 3, where we can see that, as for the 
H2 S(l)/f =2— >1 S(l)] ratio, a wide range of values 

is found (the quoted values, not corrected for extinction, 
again have an estimated uncertainty of ~20%). 

Only four sources (out of ten) are found to be 
Br 7 -dominated. These are: IRAS 06556+1623, IRAS 
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Table 3. Flux ratios observed in the sources showing molecular hydrogen emission and prevalent excitation mechanism 
tentatively inferred from the data 
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07027-7934, IRAS 10197-5750 and IRAS 17119-5926. 
Detailed morphological information exists for the latter 
three sources. IRAS 07027—7934 is known to be a very 
young PN showing a round/elliptical shape. The bipolar 
morphology is only incipient in IRAS 17119—5926 which 
otherwise would look also quite round/elliptical. The only 
source in this short list with a marked bipolar morphology 
is IRAS 10197-5750 (Sahai et al. 1999), which is seen al- 
most edge-on. This source has recently been imaged with 
HST (NICMOS) in the light of R 2 [v=1^0 S(l)] (Sahai et 
al. 2000) and the spatial distribution of the emission in this 
line is found to be very similar to the distribution of scat- 
tered light in the bright lobes, suggesting that fluorescence 
excitation is the main mechanism contributing to the over- 
all observed emission. However, a clear enhancement of 
emission is observed at two localized point-symmetric re- 
gions distant from the central star which unfortunately 
have not been covered by our slit (the slit orientation was 
90°and its bipolar axis is orientated <~ +20°), which we 
interpret as shock-excited emission probably originated by 
the interaction of the high velocity post-AGB wind with 
the slow moving old AGB shell. Interestingly, the four Br7- 
dominated sources show H 2 [v=1^0 S(l)/u=2— »1 S(l)] ra- 
tios compatible with fluorescence excitation (see Table 3). 

Two other sources in our list have also been imaged 
in the light of H 2 [v=l^0 S(l)] with HST (NICMOS). 
These are IRAS 10178-5958 (Sahai et al. 2000) and IRAS 
17150-3224 (Su et al. 2000) belonging to the group of H 2 - 
dominated sources in our sample. 

In both cases, the spatial distribution of the molec- 
ular hydrogen emission suggests a shock-excited origin. 
In IRAS 10178—5958, the molecular hydrogen emission 
comes mainly from the walls of the bipolar lobes (Garcia- 
Lario et al. 1999b), while in IRAS 17150-3224 the bulk 
of the detected emission arises in an expanding molecular 
torus surrounding the central star (Weintraub et al. 1998) 
although several clumps of molecular hydrogen emission 
are also detected in the outer regions of the bipolar lobes 
according to Su et al. (2000). A similar spatial distribution 



was also found by Sahai et al. (1998) in the case of AFGL 
2688, where the emission observed was also interpreted as 
due to shock-excitation. 

Unfortunately, none of the four other H 2 -dominated 
sources in Table 3 have been imaged in molecular hydro- 
gen with HST. Three of them (namely IRAS 14331-6435, 
IRAS 17311-4924 and IRAS 18062+2419) are sources 
not spatially resolved from ground-based observatories. 
Only IRAS 16594—4656 was found to be extended from 
the ground in optical CCD images and it was later im- 
aged with HST in the optical showing an extended, com- 
plex multi-polar morphology (Garcfa-Lario et al. 1999a; 
Hrivnak et al. 1999). All of them show large H 2 [u=1^0 
S(l)/w=2^1 S(l)] ratios which are consistent with shock 
excitation. 

The emission seems to be extended only in the case of 
IRAS 10178-5958 and IRAS 17119-5926 but we did not 
find any spatial dependence of the measured H 2 [u=1^0 
S(l)/v=2^1 S(l)] ratios along the nebulae. This con- 
firms shock-excitation as the prevalent mechanism of ex- 
citation of the molecular hydrogen in IRAS 10178—5958 
throughout its whole nebular extension. On the other 
hand, fluorescence-excitation seems to be the mechanism 
through which molecular hydrogen gets excited in the case 
of IRAS 17119-5926. For IRAS 17119-5926, 2D spectra 
of both H 2 lines were taken with different (perpendicu- 
lar) slit orientations in 1993 and 1994. Remarkably, the 
molecular hydrogen emission seems to be quite homoge- 
neous independent of the slit orientation considered and 
extends well beyond the limits of the Br7 emission. This 
suggests that most of the emitting molecular hydrogen is 
located outside the small ionized region of this very young 
PN, in the still neutral outer shell. 

In summary, when we integrated all the flux received 
in each line, we obtain high H 2 [v=1^0 S(l)/u=2->l S(l)] 
and H 2 [v=l-^0 S(l)]/Br7 ratios that are consistent with 
the prevalence of shock excitation, while low ratios in- 
dicate fluorescence as the dominant mechanism for the 
observed H 2 emission. This is supported by the different 



Garcfa-Hernandez et al.: Near— IR spectroscopy of planetary nebulae precursors 



13 



spatial distribution seen in those objects with existing H2 
NICMOS images. 

With our small sample and the rest of available data 
we find different types of sources that exhibit H2 emission 
during the post-AGB phase. We propose that this scenario 
could be summarized as follows: 

- In some cases, especially strongly bipolar sources with 
known early spectral types, there must be a contribution 
from the two excitation mechanisms in the observed emis- 
sion (e.g. IRAS 10197-5750). It is very difficult to distin- 
guish the contribution of both excitation mechanisms but 
it seems clear that the shocked H2 emission is more intense 
than the fluorescence-excited emission, leading to higher 
R 2 [v=1^0 S(l)/v=2-»l S(l)] and H 2 [w=1^0 S(l)]/Bry 
ratios. 

- The H 2 emission that comes from bipolar sources 
with later than A spectral types (e.g. AFGL 2688 and 
IRAS 17150—3224) is shock-excited because the central 
post-AGB star is not hot enough to produce fluorescence 
emission, which, so it has been demonstrated, is corre- 
lated with the evolutionary stage (spectral type) of the 
central post-AGB star. So the presence of shock-excited 
emission early in the post-AGB phase must be indicative 
of the evolutionary morphological stage of the nebulosity 
that surrounds the central post-AGB star. However, it re- 
mains unclear if the early development of marked bipolar 
morphologies during the post-AGB phase could be depen- 
dent on the mass of the central progenitor star. For higher 
masses, the morphological evolution could be faster and 
the formation of bipolar morphologies is easier; so, there- 
fore, is the detection of shock-excited H 2 emission. 

- On the other hand, the fact that H 2 emission is de- 
tected towards well evolved bipolar PNe could be indica- 
tive that the destruction of neutral material around the 
non-bipolar ones (e.g. IRAS 07027-7934 or even IRAS 
17119—3224) can roughly coincide with the photoioniza- 
tion of the circumstellar shell. This could indicate that the 
fluorescence-excited H 2 emission detected in non-bipolar 
sources is observed only during a short period of time in 
the post-AGB phase. However, in those bipolar objects 
where the wind is strongly collimated by the presence of 
an optically thick equatorial disc or torus, the molecular 
material can be shielded from the ionizing photons coming 
from the central star. 

This entire hypothesis deserves further study to under- 
stand better the physical and chemical processes involved 
during the short transition phase from the AGB phase to 
the PN stage. At present, the number of known transi- 
tion objects and the obtaining of high quality data are 
growing continuously. A detailed spatially resolved spec- 
troscopic analysis of the H 2 emission line spectrum in a 
large sample of planetary nebulae precursors is needed to 
address these issues fully. 



6. Conclusions 

With the data here presented we have extended the num- 
ber of post-AGB stars and/or proto-PNe where molecular 
hydrogen emission has been detected from four to thirteen. 

Our results confirm the previous finding by Weintraub 
et al. (1998) made on a smaller sample of sources that the 
onset of emission from molecular hydrogen takes place in 
the post-AGB phase. In many cases this occurs well after 
the generation of the bipolar structure takes place and 
before the nebular envelope gets ionized. 

In the past it has been suggested that molecular hydro- 
gen emission is a signpost for bipolarity in PNe (Kastner 
et al. 1996). We find that the detection of molecular hy- 
drogen strongly depends on the excitation mechanism. 

The detection of fluorescence-excited molecular hydro- 
gen emission in precursors of PNe is strongly correlated 
with the evolutionary stage of the central star (spectral 
type), independent of the morphology of the source con- 
sidered. The fluorescence excitation of molecular hydrogen 
emission is produced as a consequence of the absorption 
of the increasing number of UV photons escaping from the 
relatively hot, rapidly evolving central post-AGB star by 
the hydrogen molecules present in the circumstellar shell. 
As a general trend, we find that fluorescent H 2 emission 
becomes active when the central post-AGB star reaches 
a certain temperature that corresponds to an A-spectral 
type. 

On the other hand, shock-excited emission is only de- 
tected in objects with a marked bipolar morphology, some- 
times at a very early stage during the post-AGB phase, 
when the star is still very cool and, sometimes, localized 
in the waist and/or in the specific regions of the bipolar 
lobes (either at the walls, like in IRAS 10178-5958, or 
at the tips, like in IRAS 10197—5750), where a physical 
interaction exists between fast- and slow-moving material. 

Bipolarity also helps preserving the molecular hydro- 
gen from destruction, favoring its detection in well evolved 
bipolar PNe. Actually, with very few exceptions, the 
brightest H 2 emission, shock-excited, is detected toward 
the waist of bipolar PNe (Kastner et al. 1996). The high 
correlation found between bipolarity and molecular hydro- 
gen emission seems to be just the consequence of the bias 
induced by this fact. 

The diagnostic ratios H 2 [v=1^0 S(l)/i;=2->l S(l)] 
and H 2 [w=1^0 S(l)]/Br7 are shown to be very good in- 
dicators of the prevalent mechanism of excitation in proto- 
PNe. Unfortunately, the information on the relative spa- 
tial distribution of the different molecular hydrogen lines 
with respect to the hydrogen recombination lines is still 
sparse. Without this information we cannot determine 
the precise localization of shock-excited and fluorescence- 
excited regions on individual sources. As we have shown 
through this paper this is a strong limitation if we want to 
understand the physical processes which are taking place 
in this still poorly known transition phase which precedes 
the formation of PNe. This is especially important to 
interpret the emission observed in complex objects like 
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IRAS 10197—5750, where both fluorescence- and shock- 
excited regions may coexist. High spatial resolution im- 
ages (e.g. as in the past with NICMOS on board HST) 
in the three emission lines above mentioned for a statis- 
tically significant number of transition objects and young 
PNe will provide this crucial information. This would be 
an ideal way to test the scenario here proposed. 
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